home *** CD-ROM | disk | FTP | other *** search
/ QRZ! Ham Radio 8 / QRZ Ham Radio Callsign Database - Volume 8.iso / pc / files / dsp / 96000tar.z / 96000tar / 96000 / appb / b117.asm < prev    next >
Assembly Source File  |  1992-04-28  |  3KB  |  82 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.17    FIR Lattice Filter  
  8. ;N refers to the number of 'k' coefficients in the lattice filter.  Some  filters may have other coeffi-
  9. ;cients other than the 'k' coefficients but  their number may be determined from k.  
  10. ;                                          FIR LATTICE FILTER 
  11. ;       
  12. ;     COEFFICIENT AND STATE VARIABLE STORAGE 
  13. ; 
  14. ;       R0                 R4 
  15. ;                            
  16. ;                          
  17. ;   x:  S1 S2 S3 Sx     y: k1 k2 k3 
  18. ; 
  19. ;        M0=3 (mod 4)       M4=2 (mod 3) 
  20. ; 
  21. ; 
  22. ;                                  SINGLE SECTION 
  23. ; 
  24. ;                
  25. ;                    t                t'     equations: 
  26. ;                                            t'=s*k+t, t'?t 
  27. ;                           k                 s'=t*k+s 
  28. ;                         
  29. ;                        
  30. ;                       
  31. ;                      
  32. ;                     
  33. ;                             k
  34. ;                            
  35. ;               
  36. ;                   s             s' 
  37. ;                                            DSP56000 IMPLEMENTATION 
  38. ; 
  39. ;                                                              Program        ICycles
  40. ;                                                              Words 
  41. ; move      #state,r0     ;point to state variable storage 
  42. ; move      #N,m0         ;N=number of k coefficients 
  43. ; move      #k,r4         ;point to k coefficients 
  44. ; move      #N-1,m4       ;mod for k's 
  45. ; 
  46. ; movep     y:datin,b     ;get input 
  47. ; 
  48. ; move      b,x:(r0)+  y:(r4)+,y0  ;save 1st state, get k    1       1 
  49. ; do        #N,_elat               ;do each section          2       3 
  50. ; move          x:(r0),a   b,y1    ;get s, copy t for mul    1       1 
  51. ; macr y1,y0,a             a,y0    ;t*k+s, copy s            1       1 
  52. ; macr x0,y0,b  a,x:(r0)+  y:(r4)+,y0 ;s*k+t, sv st, nxt k   1       1 
  53. ;_elat 
  54. ; move       x:(r0)-,x0 y:(r4)-,y0 ;adj r0,r4 w/dummy loads  1       1 
  55. ; 
  56. ; movep      b,y:datout    ;output sample                  -----   ----- 
  57. ;                                                   Totals:  7      3N+5 
  58. ; 
  59. ; 
  60. ;                                            DSP96002 IMPLEMENTATION 
  61. ; 
  62. ;                                                            Program    ICyc
  63. ;                                                            Words 
  64.  move     #state,r0     ;point to state variable storage 
  65.  move     #N,m0         ;N=number of k coefficients 
  66.  move     #k,r4         ;point to k coefficients 
  67.  move     #N-1,m4       ;mod for k's 
  68.  
  69.  move      y:datin,d5.s  ;get input 
  70.  
  71.  move                      d5.s,x:(r0)+ y:(r4)+,d4.s ;sv s,get k  1    1 
  72.  do        #N,_elat                                  ;do filter   2    3 
  73.  fmpy d5,d4,d3             x:(r0),d0.s               ;t*k, get s  1    1 
  74.  fmpy d0,d4,d1 fadd.s d3,d0                          ;s*k,t*k+s   1    1 
  75.                fadd.s d1,d5 d0.s,x:(r0)+ y:(r4)+,d4.s ;s*k+t; s,k 1    1 
  76. _elat 
  77.  move         x:(r0)-,d0.s y:(r4)-,d7.s ;adj r0,r4 w/dummy loads  1    1 
  78.  
  79.  movep        d5,y:datout            ;output sample 
  80. ;                                                                ---  --- 
  81. ;                                                        Totals:  7  3N+5 
  82.